Existence of a Density Functional for an Intrinsic State 



B. G. Giraud 
bertrand.giraud@cea.fr, Institut de Physique Theorique, 
DSM, CE Saclay, F-91191 Gtf/Yvette, France 

B. K. Jennings 
jennings@tnumf.ca, TRIUMF, Vancouver BC, V6T2A3, Canada 

B. R. Barrett 
bharrett@physics.arizona.edu, Department of Physics, 
University of Arizona, Tucson, AZ 85721, USA 
(Dated: August 19, 2008) 

A generalization of the Hohenberg-Kohn theorem for finite systems proves the existence of a 
density functional (DF) for a symmetry violating intrinsic state, out of which a physical state with 
good quantum numbers can be projected. 



I. INTRODUCTION 

Density functional ^ theory (DFT) was initially defined for ground states. These have good quantum numbers. 
Every nuclear physicist knows that, for instance, the ground state of ^"^Ne is a 0'*' and that its density is, thus, isotropic, 
not an ellipsoid. Every molecular physicist knows that, for instance, the ground state of the ammonia molecule is a 
good parity state, not just the pyramid described by the Born-Oppenheimer approximation. In particular, the nuclear 
DF must generate .spherical solutions for the some thousand 0"*" nuclear ground states, whether nuclei are intrinsically 
deformed or not. The same need for isotropic solutions extends to the non-local generalization of the DFT - [1]. 
But the theory of rotational bands and/or parity vibrations, whether in nuclear, atomic or molecular physics, most 
often relates ground states to wave packets, often named "intrinsic states" , which are symmetry breaking, namely do 
not transform in an irreducible representation (irrep) of the symmetry group S of the Hamiltonian. Therefore, one 
may raise the question of DFT for intrinsic states rather than eigenstates. 

Given the physical Hamiltonian H with its symmetry group S, calculations providing a "non S-irrep" state as a 
solution for a minimum energy cannot be labelled as the result of "the" DF. Such a state, labelled intrinsic, is actually 
just a convenient wave packet, to be subsequently projected onto good quantum numbers to account for physical 
levels. Such intrinsic calculations should rather exhibit a special Hamiltonian, which might be called an intrinsic 
Hamiltonian, distinct from the physical one, if such calculations are to be legitimized. Or they should be interpreted 
as one variety of the Hartree-Fock, Hartree-Bogoliubov, etc. variational approaches. This is implicit or even explicit 
in calculations with an energy density functional, implying non- localities, see for instance [J] - |7[. Energy density and 
particle density are different concepts. 

It turns out that the particle density which has been used for the foundation of DF theory mainly concerns eigen- 
states of the physical Hamiltonian, in principle at least, while the energy density, used for Skyrme force calculations 
in nuclear physics for instance, mainly provides intrinsic states. This note presents a particle DF theory for intrinsic 
states, not for eigenstates of the Hamiltonian. We show how the physical Hamiltonian can be reconciled with the 
proper definition of a DF for an intrinsic state and how the resulting intrinsic state can be accepted as a useful wave 
packet, out of which states with good quantum numbers can be projected. In particular, in the case of molecules, 
our approach will consider both the electrons and the nuclei. Our intrinsic state can take into account both kinds of 
degrees of freedom. Section II describes a functional out of which a variational principle derives for an intrinsic state, 
and out of which a DF for the intrinsic density is obtained. Section III gives an example of variational equations to 
be solved in practice. Section IV rewrites the formalism into a slightly simpler form. Section V contains a discussion 
of our result and suggests an ansatz for intrinsic Hamiltonians. 



II. BASIC FORMALISM 



For a first argument, dealing with one kind of identical particles only, let H be their physical Hamiltonian and 
0, (010) = 1, be a trial wave packet, most often not transforming under an irrep of the symmetry group S of H. 
For instance, for fermions, may be an arbitrary Slater determinant, but we let be also a more general wave 
function, including some amount of correlations. States ip oc P0 with good quantum numbers can then be projected 
out of by a projector P, a fixed operator. In the following, we shall systematically use the properties, = P and 
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[P, H] = 0. It may happen that {(j)\P\(j>) vanishes, but such cases usually make a domain of zero measure in the usual 
variational domains, where (f> evolves. In any case, since H is an operator bounded from below, the functional of 0, 
{(j)\PH\(j>) / {(j)\P\(j)) , is bounded from below. Embed now the system in an external, local field, U = The 
local, real potential u is taken bounded from below, but is otherwise arbitrary. In particular, it may usually have 
none of the symmetries of H. Then, given u, the following functional of (p, 



\PH[ 



(1) 



is bounded from below. To find the lowest energy with the quantum numbers specified by P one can use a constrained 
search [S], in which one first considers only states that show a given density profile T(r), then one lets r vary. 
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The process goes in two steps, namely, i) a minimization within a given particle density profile, T(r) = ((/)|cJ!cr|(/>), for 
N particles, with cl and Cr the usual creation and annihilation operators at position r, then, ii) a minimization with 
respect to the profile. The inner minimization clearly defines a DF, F[t\ = Inf^^T- {{4)\PH\4)) / {(j)\P\(j))) . 

Actually, it is more general Q - [l^ to use many-body density matrices B in A^-body space, meaning mixed as well 
as pure states, and yielding a density T(r) in one-body space. 
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but we shall use wave-functions in the following, namely B — \(f>){(j)\, for obvious pedagogical reasons. We shall assume 
that this Inf^ actually defines an absolute minimum, Min^, reached at some solution $ of the corresponding variational 
principle. Moreover, we shall assume, temporarily at least, that the solution $ is unique. Uniqueness is not obvious, 
however, if only because many (/)'s can give the same P\(t>) , and, when u vanishes, this variational principle, Eq. 
reduces to the well-known "variation after projection" [lj| method for Hartree-Fock calculations for instance. 

Anyhow, r and u are clearly conjugate in a functional Legendre transform, with SF/St = —u. Finally, if p{r) 
denotes the profile of $ when u —>■ 0, then the lowest energy with good quantum numbers is nothing but F[p\. 
The minimization, with respect to r, of the functional, F[t], provides simultaneously the density of the intrinsic 
(unprojected!) state and the projected energy. Notice, incidentally, that F[t] depends on the choice of the variational 
set of trial functions where the "inner minimization" is performed. Furthermore, it obviously depends on P. 

A more general argument is possible, with more than one kind of identical particles. Trial states can be, for instance, 
products of determinants, one for each kind of fermions, and permanents, one for each kind of bosons. Consider for 
instance the ammonia molecule, with i) its active electrons, ii) its three protons and iii) its nitrogen ion. It is trivial 
to include a center-of-mass trap into H to factorize into a spherical wave packet the center of mass motion of this 
self-bound system and avoid translational degeneracy problems. The complete Hamiltonian trial states </> and 
density operators B depend on and describe simultaneously the electron, proton and nitrogen ion coordinates and 
momenta. A DF in just the electronic density space, however, results from the definition. 



T[t] = Inf£ 
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where P projects good quantum numbers for the whole system and traces are taken over all degrees of freedom, while t 
is set as only an electronic density. Interactions between heavy degrees of freedom, between heavy and electronic ones, 
and between electrons, are taken into account by the trace in the numerator. No Born-Oppenheimer approximation 
is needed for this "global" definition. For the sake of simplicity, however, we return in the following to the case of one 
kind of particles only. Most considerations which follow have obvious generalization for multicomponent systems. 



III. VARIATIONAL EQUATIONS 



Let d(j) be an infinitesimal variation of the trial function in its allowed domain. Then, at first order, one obtains, 

'(b\PH\(h) {5<h\P\(h) {6\PH\6) 
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If one defines the "gradient operator" , 

G = + r/ - P^'t>\PH\<j^) 

(^|P|0) ^ mp\<p)? ' ^ ^ 

then, obviously, 51F — {5(t)\Q\4)) + ((/!)|f/|(S0). Note, incidentally, that Q is Hermitian. 

At the minimum position the variation 5J- vanishes for any 5<j}. Replace 5(j) by iScf) to see that the difference, 
— {5(j)\Q\^) + ($|tj|(50), vanishes as well. Then, trivially, at <i>, both ((5(/>|tJ|$) and {i^\Q\5(t>) vanish simultaneously, 

(<S0|e|$) = ($15150) = 0, V(50. (7) 

In the special case of Slater determinants, let \ph) = cj,Ch\^) denote any particle-hole state built upon |$) as the 
"reference vacuum" for quasi-particles. Here and c are the familiar fermionic creation and annihilation operators, 
respectively. Then the particle-hole matrix elements of Q vanish, 

{ph\ g 1$) = 0, Vph. (8) 

As long as a solution of this stationarity condition, Eq. ([H]), is not reached, the matrix elements, {ph\Q\(p), define 
the direction of the gradient of J- in the hyperplane tangent to the manifold of Slater determinants. A gradient 
descent algorithm, \d(p) — ~vJ2ph \P^) (p^\G\'P)j where 77 is a small step parameter, then leads to the solution. Notice, 
however, that the ph representation is covariant with 0. The ph basis has to be recalculated at each step. Being state 
dependent, Q must also be recalculated at each step. 



IV. SIMILAR THEORY, WITH A LAGRANGE MULTIPLIER 

The slightly complicated gradient operator, Eq. ([6]), leads to a variational condition, Eq. ([7]), which combines 
the matrix elements of three operators, namely PH, P and U. Define the number A = |$)/($|P|$) as a yet 

unknown Lagrange multiplier; it can be considered as an arbitrary parameter and shall be adjusted self-consistcntly 
later, when <I> is reached. Then Eq. ([7]) also reads, 

{6(j)\ {PH - XP + {<^>\P\<i>)U) \<^>) = 0. (9) 

If (j) were completely unrestricted, namely if dcj) were completely general, this equation, Eq. ([5]), would mean that $ is 
an eigenstate of the operator Q. Since intrinsic states are understood to belong to restricted sets of states, the result 
$ is only an approximate eigenstate of Q. 

To avoid the cumbersome coefficient, ($|P|$}, which multiplies U, it is convenient to define the auxiliary operator, 

H = PH -XP + W, (10) 

where W = J2i w{ri) is, like U, an arbitrary, local, real, external field, bounded from below. It is obvious that Q and 
Ti define a common solution $ if u and w are suitably proportional to each other, w = ($|P|$) u. In the following, 
however, we set H. ab initio. It is an operator bounded from below. We are interested in its "almost ground state" 
S and assume that this state is unique. A connection between a solution S in this section and a solution $ in the 
previous section can easily be tested later. 

Define again a constrained search for the lowest energy, 

H|0)=Inf,(Inf^^,(0|7^|0})-Inf, (^PA[T]-t- J drT{r)w{r)Y (11) 

where the A-dependent DP, Fx, is defined as, 

Pa[t] =Inf^^, (01 (Pi/- AP) !</)). (12) 

It is again convenient, for pedagogy at least, to assume that this Inf^ induces an absolute minimum, reached at a 
position S in the variational space. The same assumption states that, given A, the absolute minimum of Fx[t] is 

Fx[a]^{E\{PH-XP)\E), (13) 

where a is the density of S. Let £ denote this energy, £{X) = (S|(Pii — AP)|S). A simple manipulation then gives, 

| = -(Him- (») 

A Legendre transform, using A and (^|P|,^) as conjugate variables, is thus available to return the matrix element 
(S|Pii|S) as a function of the matrix element (S|P|S). Then one has just to locate the minimum of their ratio. 

Note again that the theory depends on the variational space where 4> evolves. But, in any case, one obtains 
simultaneously the density of S, the best intrinsic state, and the energy of its projected state P| S). 
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V. SUMMARY, DISCUSSION AND CONCLUSION 

If only because of the need for spin densities (isj in the description of polarizable systems, the problem of symmetry 
conservation, or restoration, in DF theory has already received much attention in atomic and molecular physics [l6[ 
-_[l9|. It has been revisited here, in the spirit of the projected Hartree-Fock method with variation after projection 
[14| : a variational principle for the density of an intrinsic state, without symmetry, optimizes the energy of a state 
with good quantum numbers. The idea was already introduced in the context of particle number projection [20j . We 
have shown in Sees. II and IV that our approach allows generalizations of the Hohenberg-Kohn existence theorem. 

It can be stressed that the present approach is concerned with the density of an intrinsic state, not that of an 
eigenstate. This is a major difference with all the other DF theories that we are aware of. Note, in particular, how 
our functional differs from a functional of a symmetrized [l6l | density . 

We showed in Sec. IV that a way to define the intrinsic Hamiltonian amounts to a linear combination, 7i = 
—A P + PH, of the projector P on the desired quantum numbers, and the laboratory Hamiltonian multiplied by that 
same P. Here, a subtle question must be raised, that of the nature of the intrinsic state. The more flexible the trial 
functions for this state, the better the projected state and the lower the projected energy. However, full flexibility 
contradicts simplicity, and, moreover, uniqueness of the intrinsic state; many different packets can give the same 
P|(/)). Symmetry projection brings correlations which, therefore, must be absent from the intrinsic state. This is why 
variational domains for intrinsic states must necessarily be much narrower than the full Hilbert space. 

In practice, fortunately, intrinsic states are confined to non-linear, curved [2ll | manifolds, such as coherent states, 
Slater determinants, etc., which do not make linear subspaces. The intrinsic state, therefore, is not an exact eigenstate 
of 7i. It just minimizes a related quantity, the projected energy. It must be concluded that DF theory for an intrinsic 
state necessarily depends on two factors, namely, i) obviously the quantum numbers to be projected out, but also ii) 
the variational space retained for this intrinsic state. 
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